A Vishniac type contribution to the polarisation of the CMBR? 
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Abstract 



pLn ', Radiation which has a quadrupole component of anisotropy, can get po- 

Q"^ ' larized by Thomson scattering from charged particles. In the cosmological 

context, the microwave background photons develop significant quadrupole 
^ . anisotropy as they free stream away from the the epoch of standard recombi- 

0^ , nation. Reionization in the post recombination era can provide free electrons 

to Thomson scatter the incident anisotropic CMBR photons. We compute 
' the resulting polarisation anisotropy on small (arc-minute) angular scales. 

' We look for significant non-linear contributions, as in the case of Vishniac 

0> ■ effect in temperature anisotropy, due to the coupling of small-scale electron 

density fluctuations, at the new last scattering surface, and the temperature 
{pH| quadrupole. We show that, while, in cold dark matter type models, this does 

O ' not lead to very significant signals (~ 0.02 — 0.04/ii('), a larger small angular 

' scale polarization anisotropy, (~ 0.1 — 0.5/x/C), can result in isocurvature type 

■ models. 



I. INTRODUCTION 

The cosmic microwave background radiation has a wealth of information about the pa- 
rameters which govern the dynamics and the physical processes in the universe. Three 
important aspects of the CMBR in which the information about the universe is encoded are 
its temperature anisotropy, spectral distortion and anisotropy in polarization. A study of 
these aspects will allow one to constrain the evolution of both the background universe and 
the large scale structures. 

Polarization of the CMBR arises from the Thomson scattering of radiation from free 
electrons. An analysis of polarization properties of the CMBR can reveal valuable informa- 
tion about the ionization history of the universe. Much of the current work on polarization 
anisotropy is in the context of the linear perturbation theory. The expectation generally 
is that higher order terms are much smaller. However, for the temperature anisotropy, on 
very small angular scales, of order of arcminutes, it was shown by Vishniac that nonlinear 
effects can also make significant contributions These arise through mode coupling of 
the electron density perturbations on small scales with source terms which vary over larger 
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scales. The Vishniac effect is especially important in models where there is significant early 
re-ionisation ( |^). 

Zaldarriaga 0] studied the effects of such early reionisation on first order polarisation 
anisotropy, in a semi-analytical fashion. It turns out that the CMBR can develop a significant 
quadrupole, by the epoch of re-ionisation, due to the free streaming of the monopole at 
recombination. The Thomson scattering of this quadrupole, off the electrons at the re- 
ionised epcoh, can lead to additional polarisation signals. In this paper we wish to follow 
suit and examine whether this quadrupole, coupling to fluctuations in the electron density, 
at the new last scattering surface, can also lead to significant, Vishniac type second order 
efects. And result in a polarization anisotropy of the CMBR at small angular scales. 

In the next section, we give the basic equations and their formal solution. Section 
III presents an analytical estimate of the Vishniac type contribution to the small scale 
polarization anisotropy, in re-ionised models. Section IV gives numerical values for some 
illustrative models of structure formation. We summarise our conclusions in Section V. 



II. BASIC EQUATIONS AND THEIR FORMAL SOLUTION 

The equations governing the evolution of polarization perturbation Ap(x, 7, r) and tem- 
perature perturbation At(x, 7,r) = AT/T, for scalar modes can be derived from the mo- 
ments of the Boltzman equation for photons. In the conformal Newtonian gauge, they are 
given by 

Ap + jAAp = neara(r)(-Ap + ^[1 - P2{fi)]Il) (2.1) 

At + l^^iAT = (/> - lid^ip + n,aTa{T){-AT + + liVi - ^P2(/x)n) (2.2) 

Here x is the comoving co-ordinate, r is conformal time. He the electron density, v the 
fluid velocity, 7 is the direction of photon propogation, and ip the conformal Newtonian 
potentials, and a dot represents derivative with respect to conformal time. We have also 
defined 

n(x, r) = -At2(x, r) - Ap2(x, r) + Apo(x, r). (2.3) 

with Ato, Apo the monoplole, and At2, Ap2 the quadrupole temperature and polarisation 
perturbations, respectively. We define these angular moments by 

Ap,t(x, 7, r) = ^(2/ + l)PK/i)A«,TKx, r); Api^ri = j ^Piifi)Ap,Tifi) (2.4) 

[Here, we make the usual assumption that after Fourier transformation, the evolution equa- 
tions for temperature and polarisation perturbations, depend on 7, only in the combination 
7.k, where k is the wavevector; so we define /i = j.h/k, where k = |k| (cf. ref. p|)]. 

In order to treat inhomogenieties in the electron density, we take ne(x, r) = ne(r)[l + 
(5e(x, r)], where 6e is the fractional perturbation of electron density about the space averaged 
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mean. In general we will have Se << 1- We may note at this point that spatial perturbations 
in the number density of the electrons is precisely the feature that gives rise to Vishniac 
effect in second order temperature perturbations. We will investigate a similar effect for 
polarization perturbations. 



It will be convenient to express equation ^TT| in terms of the fourier modes as follows, 
Ap + ikfiAp = neO"Ta 



^(1 - P2(^))n(k, r) - ^(1 - P2(/x))5(k, r) - Ap 



(2.5) 



We have retained the same symbols for the fourier transformed quantities and defined the 
mode coupling source term S{k, r) by 

Sik, r) = -j ^5e(k - p, r) [n(p, r) - Ap(p, r)] (2.6) 

The formal solution of Eq. ( p.5|) is given by 

Ap(k, r)= rUl- P2(/i)) [n(k, r') - 5(k, r')] g{T, ry'^^^'-^Ur' (2.7) 
Jo 2 

where the g{T, r') called the visibility function is given by 

g(r,T') = ne(r')(TTa(r')e-/>^("")"^"("")''"" (2.8) 

In the above equations, the value of the polarization perturbation at the epoch r is deter- 
mined by the entire history from r' = to r' = r. The visibility function g{T, r') determines 
the probability that a photon last scatterred at the epoch r' reaches us at the epoch r. 
The exact form of the visibility function is determined by the ionization history of the uni- 
verse. Operationally, the role of the visibility function is to give different weightages for the 
integrand for different epochs. 

In this paper we consider a model in which the universe underwent a phase of standard 
recombination and got reionized completely at a later epoch, t^. We will only be concen- 
trating on the second order polarization perturbations arising from the ^(k, r) term in Eq. 
( p.7| ). The S'(k, r) contribution to the RHS of Eq. ( |2.7|) involves a convolution in the fourier 
space, which couples the first order temperature (polarisation) perturbations with the first 
order perturbations in the electron density. A very similar situation exists in the case of 
Vishniac effect in second order temperature perturbations. So we expect a Vishniac type 
effect in second order polarization perturbations as well. 

Further, the coupling of Ap, Ap,, and Ap,^ with Se in S, are likely to be much smaller than 
the coupling of 6e and At2- This is because firstly the temperature perturabtions generically 
dominate the polarisation. Also the quadrupole temperature anisotropy, Ar2(k, r), will grow 
to a larger value, between the epochs of recombination and reionisation, due to free streaming 
of the monopole at recombination (cf. [Q). We therefore retain in S only the Aj'2 term and 
neglect the other terms. We then have 

We will adopt this approximate form for the mode coupling term in what follows. 
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III. VISHNIAC TYPE CONTRIBUTION IN REIONIZED MODELS: ANALYTIC 

ESTIMATE 



The absence of " Gunn- Peterson" dips in the spectra of distant quasars indicates that the 
universe was probably reionised at some redshift z = z^, > 5. (0). The value of is not 
known observationally, while different theoretical models have different predictions for this 
redshift. In the model which we consider, the universe underwent standard recombination 
at Tr, and was reionized completely at a later epoch t^. In this case as shown in Ref. Q], the 
visibility function has two peaks, one around and another peak around r^. We wish to 
consider here Vishniac type, second order contribution to Ap. This comes dominantly from 
the value of r' around the latter peak. It is then convenient to separate the integral over r' 
in Eq. (|2.7|) in two parts: < r' < r^, and r^, < r' < Tq, where Tq is the present conformal 
time. We write Ap = Ap + Ap, with 

A^k, r) = r PM) [n(k, r') - 5(k, r')] g{r, ry'^^^'-^^dr' (3.1) 

A^(k, r) = 1(1 - P2(/i)) [n(k, r') - 5(k, t')\ g{T, rY^^^^' dr' (3.2) 

The first contribution in ( p.l| ) is simply Ap = exp (— K*)Ap'^, Here Ap^ is the polari- 
sation that would be measured if there was no reionisation and is the optical depth to 
Thomson scattering between now and recombination. This contribution is reduced by the 
the fact that only a fraction exp (— k^,) of the photons that arrive at the observer come di- 
rectly from recombination, without further scattering. (Also the second order contributions 
from the S term are much smaller than the first order term because the electron density 
fluctuations at recombination 5e{Tr) ~ 10~^ — 10"'^ << 1, for the relevant k values). 

In order to calculate the second contribution, one has to determine the form of the 
visibility function after the standard recombination epoch, that is g{TQ, r') for t' > t^. Using 
the exact form for g{TQ, r'), to solve for the Ap is not analytically tractable. So we resort to 
an approximation for g{To, t') in this work, which, while preserving its main features, also 
allows analytical results to be derived. We will return to a full numerical treatment of the 
problem elsewhere. In particular, we choose the form of the visibility function after standard 
recombination, to be a truncated exponential, given by 

g(ro,T')=N-e-^9iT'-n) (3.3) 
a 

Here the Heavyside 6{x) function, is zero for x < and 1 for x > 0. It takes account of the 
fact that before reionisation, rig = 0. Further, is a normalisation consant and a gives the 
spread of the exponential. By appropriately chosing a, we can set the width of the reionised 
last scattering surface. Also note that g{To,T') has the interpretation of probability; so its 
integral over r' from r' = to r' = tq should be normalized to unity. This determines 
the normalisation factor A^. For a sufficiently early epoch of reionisation, we generally have 
(tq — T^)/(y2 » 1- In this case, the condition that the integral of g{TQ, r') over r' should be 
unity implies + e"*^* = 1, or = 1 — e""*. So A^ measures the probability of at least one 
scattering between tq and r^,, due to the reionisation. Another feature to note is that in our 
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approximation, tq does not appear at all. This is because, for the models we will consider, 
the major contribution to the scattering optical depth comes from epochs much before Tq. 

In the equation ( |3.2| ) for Ap, the first order contribution to the polarisation due to a 
reionised universe has already been considered in detail by Zaldarriaga Q. So, here, we 
concentrate on purely the second order Vishniac type effect, due to 5*, which we call Ap. 
This can be written as 

A^ = -^[1 - P2(/^)] r exp^'=^(^'--) giro, r')S{k, r')dr' (3.4) 



where retaining only the At2 term the mode coupling term is given by Eq. (|2.9| ). In 
evaluating the r' integral in Eq. (|3.4|) we assume the visibility function to be given by the 
tuncated exponential form of Eq. ( p.3|) . 

Let us look at the mode coupling term ^(k, r), given by Eq. (|2.9|) , in a little more detail. 
This term involves a coupling of the quadrupole temperature perturbation at r > r^,, and 
the electron density perturbation at the same epoch. Note that the temperature quadrupole 
at late times, can have a significant contribution due to the free-streaming of the monopole. 
For example, in a flat universe, at large enough wavelengths, the first order quadrupole 
temperature perturbation, is related to the temperature perturbations at recombination by 



At2(P, Te) = [Ato(P, Tr) + 1p{p, r^)]j2b(r - T,.)] (3.5) 

where j2 is the second order spherical Bessel function, and p = |p|. (Here we have assumed 
that p is small enough that the doppler velocity term makes little contribution to the fre- 
streamed qudrapole cf. 0). The mode coupling term can then be written as 

^^^5e(k - p, r)[ATo(p, Tr) + ^(p, Tr)]j2[piT - Tr)]. (3.6) 

The spherical Bessel function, j2{0y is well approximated by a gaussian peaked at ^ ~ 3.345 
and with a spread of ~ 1.4. So j2[p{T — Tr)] is peaked at wavenumbers around p = Po ^ 
3.345/(r — Tr). Note that for the small scale polarisation anisotropy which we wish to 
calculate, k = |k| >> po; in general we will have k ~ {10Mpc)~^ to {lMpc)~^, where 
as Po ~ (300Mpc)-^ to (1000Mpc)-\ for r > r, (cf. Ref. @] and see below). So in the 
mode coupling integral, for a fixed k >> po, the electron density perturbation 5e((k — p), r) 
varies negligibly with p, in the range of p for which the the bessel funtion makes significant 
contribution. So we can evaluate 6e at p = po and pull it out of the p integral. Also since 
k >> Po, one can approximate (k — po) ~ k. The mode coupling ingral, for large k » po 
then simplifies to the uncoupled form, 

S{k, t) = 4(k, r) J ^At2(p, r) = 5e(k, t)Q2{t) (3.7) 



where we have used Eq. ( |3.5| ) to rewrite the resulting expressions in terms of A^a again, and 
defined for later convenience, 



Q2(r)= /7^A^2(p,r) (3.8) 
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We now evaluate Ap, using Eq. (p.4|), and Eq. ( p.3|) . Let us assume that the re-ionisation 
epoch is late enough, that the electron density perturbation trace the perturbations in 
total matter density. Then, in a fiat universe, the electron density perturbations grow as 
a(r) oc r^. For an open model or one involving a cosmological constant A, the growth law 
is given as 

(3.) 

Here, tq is present conformal time and the functions / takes into account the reduced growth 
in the open or fiat A models at late times, compared to a fiat, dark matter dominated model. 
In a fiat universe Q{z) = Qq = 1, /(I) = 1 and we recover the 6e oc growth law. The 



expressions for /, and il{z) for A dominated and the open models can be found in ref. ||T5 
and fl^, respectively (in these papers, our f{fl) is called g{fl))- For z » 1, or t/tq << 1 
and f{fl{z)) 1. So the approximation 6e{k,T) = {T/To)^6e{k,To)f^^{Qo), works very well 
for these other models as well, at sufficiently early times. 

Further, as displayed explicitely in ref. [0, the power in the CMB monopole, per unit 
logrithmic inteval of p space, p^[At(,(p, r^) + ip{p,Tr)] is roughly constant on scales p < 
{100Mpc)~^ . (This reflects the fact that perturbations on scales larger than the Hubble 
radius at recombination have not evolved, and are laid out with constant power, for a scale 
invariant initial power spectrum). Recall that the presence of the j2 term, in the integral 
over p, picks out dominantly contributions from p ~ po < (300Mpc)~^, at any time r > r^,. 
Now in any realisation, one does expect some variation of the monopole, as p is varied. 
Nevertheless, because of the constancy of monopole power with p, for p ^ po, one expects 
the integral term Q2{t) in Eq. ( p.7|) , to vary much slower with r, than the electron density 
perturbation 5e(k, r), or the visibility function. This will especially be so, if the visibility 
function is sufficiently peaked around the reionisation redshift (for a small enough). So 
when evaluating Ap, we will assume that Q2{t) can be evaluated at some effective Tg ~ r*, 
and pulled out of the integral over conformal time r'. 

The remaining integral, which can be done analytically, then gives 

= -^[1 - P2(/i)] exp-M-^-) J^Ii^5,(k, r.)Q2(re), (3.10) 
2 (1 — zkfia) 

where we have taken (tq — T^)/cr » 1, as before, and defined the factor F{k,fi) 

(We have also assumed r^,/ro is small enough that Q{t^) ~ 1). 

We are now in a position to compute the Vishniac type contribution to the polarisation 
power spectrum. We define this simply by analogy to the temperature power spectrum (cf. 
0). Recall that, for the temperature perturbations, one first expands in spherical harmonics, 
with 

AT 

— = At(x, 7, To) = aimYimi'j). (3.12) 

Im 
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Then the mean square temperature perturbation is 



<i^f>-Y.'-^C„^jQr{kfi (3.13) 



where 



and 



/rv (Ik 
<\ATi{k,To)\^ > (3.14) 

QT{k) = ^\ r < \At\' > rf/i. (3.15) 

(Note that, as before, we have taken the normahsation volume, over which periodic boundary 
conditions are assumed to be = 1). Qxik) gives the power in temperature perturbations 
in any logrithmic inteval of k. 

We define therefore, a corresponding Vishniac type contribution to the polarisation power 
spectrum Qp{k), given by (cf. p), 

QUk) = < l^lf > df.. (3.16) 



We calculate Qp below. For this, first take the ensemble average of |App. We get 



1 7^21^12 

<|Alir>=j|l-A(,)|=I^^J^P.ftr. 



1 /AT\ 



2 



5 V T y, 



(3.17) 



In the above we have assumed that the 6^ and are uncorrelated with each other, defined 
the pow 
defined 



the power spectrum of electron density fluctuation as, P(.{k,T^) =< |(5e(k, r^.)|^ > and also 



/AT\2 , , 5CT2(re) f d^P , . . .,2 . ^ 

Here, (AT/T)Q(re) is the quadrupole temperature anisotropy as seen by an observer at the 
conformal time Tg. 

1/ I 2 m 

Now turn to the integral of < |Ap| > over /i. The factor, 1 — P2 can be expressed as a 
sum of even powers in /x. 

[l-P2]' = lil + ^^'-2^^') (3.19) 

We are generally intersted in the behaviour of the power spectrum for large values of k, 
or small angular scales. In this case the dominant contribution to the integral over fi, in 
determining Qp, will come from the vicinity of /i = 0. It then suffices to retain only the first 
term in the above expression for (1 — P2)^- The integral over can be done analytically to 
give the remarkably simple expression 
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Qp{k,To) = -r77^G{—) 
IdU 



ka 



/AT\ 



2 



(^j (r.). (3.20) 



In doing the integral we have assumed that A;cr >> 1 is large enough that tan ^{ka) ^ 7r/2, 
and defined the factor G{y) = 1 + 2y + Sy"^ + 5y^/2 + 35?/'^/32, where y = a /t^. Further, 

is the power per until logrithmic interval in k space, of the electron density perturbations, 
at the epoch t = t^. 

We see that the contribution to polarisation anisotropy, due to the second order Vishniac 
type effect, for re- ionised models, is basically proportional to the product of the tempera- 
ture quadrupole, and the power in electron density perturbations, at last scattering. For 
small angular scales, or large k, Qp is suppressed because of the finite thickness of the last 
scattering surface (a), by a factor ka. We note that this suppression is much milder than es- 
timated in Ref [§], essentially because in that paper, the first order temperature quadrupole 
contribution due to free-streaming of the monopole at recombination, was not included. The 
power spectrum of electron density perturbations ofcourse depends on the model for struc- 
ture formation. Also the parameters a,T^,,N depend on the re-ionisation history. However 
the power in the temperature quadrupole at ~ is likely to be of order the observed 
quadrupole; for a large enough and if it arises due to the free-streaming of the monopole 
at recombination. This is again because of the slow variation of Q2{t) mentioned previously. 
We now use Eq. ( p.20|) to make numerical estimates of the polarisation due to Vishniac type 



effect in CDM and other models of structute formation. 



IV. NUMERICAL ESTIMATES IN DIFFERENT MODELS 

A. CDM and variants 

Consider first the case of a standard CDM model (SCDM), with matter density equal 
to critical density [Qq = 1), a baryonic contribution fit = 0.05, and a Hubble constant 
h = [Ho/lOOkms'^Mpc"^) = 0.5. The optical depth to Thomson scattering in a fully 
ionised, matter dominated, fiat universe is given by k{z) = OMlSQBhlil + zfl'^ - 1]. In 
general, for a universe with matter density fio) and assuming » 1, an optical depth 
is obtained at a redshift 

z.-97.1.^/3r]f(i^)-/l (4.1) 

So to have = 1, in standard CDM, we need ~ 97.1, and = 0.2 needs z^ ^ 33.2. 
Also in a fiat matter dominated universe, the conformal time is related to the redshift by 
r = ro/(l + ^)^/^ where Tq = 2H^^ = GOOOh-'^Mpc. (Note we adopt units with c=l). So 
given z^, this fixes r*; for = 1, we have r* ~ 605.8h~^Mpc, while for k^, = 0.2, we have 
^ 1025.8h-^Mpc. 

In order to estimate the parameter a, in the model visibility function ( |3.3| ), we pro- 
ceed as follows. Let = + o". From Eq. (|3.3|) , at epochs after re-ionisation, we 
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have g{jQ,T^) / gijQjTrn) = e. For the exact visibihty function, the same ratio is given by 
9iTo,n)/g{To,Tm) = [a^(rm)/a2(rj] exp(K^ - kJ, where = K(r„). Equating these 
two expressions gives an estimate of a. In particular, using k oc (1 + zY^'^, vahd for 
large z » 1 and a oc oc (1 + z)~^, we have the implicit equation for Tm/r^, 
1 = Aln{Tm/T^) — — (T*/Tm)^]- For = 1, this gives a ~ 0.54r* = 327h^^Mpc, 
while for = 0.2, one gets a = 0.32r* = 328.3h~^Mpc. 

It remains to fix Ag and the temperature quadrupole. We take the Ae(A;,ro) = A{k), 
where the matter power spectrum 

A^(fc) = ^=(^)'4(^)m (4.2) 
with the transfer function T{k) in the form given by Bardeen et al. |^ 

= lni^±2M4l 

2.34g [1 + 3.89g + (16.1g)2 + (5.64g)3 + {Q.llqY]^'^ 

where q = k/ {hV) (cf. ref. fill). The parameter F is referred to as the shape parameter, 
is given in ref. |Tl| (eq. D28 and E12). It is of order 0.48 for standard CDM. The four-year 
COBE normalisation gives Suik = Hq) = 1.94 x 10^^ (cf. ref. fl^), for a scale invariant 
initial power spectrum (with n = 1). For such initial conditions, the COBE data also give 
the {AT/T)q{to) = 18 ± l.GfiK 0. The value of n > lOOOMpc, that we generally obtain, 
is likely to be large enough so that (AT/T)Q(re) ~ {AT/T)q{tq), with reasonable accuracy. 
So we will scale the quadrupole at Te with the present day observed value. From the above 
considerations, and normalising all quantities to a value of k = kc = lhMpc~^, we get 



= 1 (4.4) 



In case = 0.2, one has to replace the numerical value in ( [4.4| ) by yQ^{kc) ~ 1.5 x 
lO^^^K. Lower values of k^, < 1 may already be implied by the observed tentative rise in 
the CMB anisotropy on degree scales. Note that decreasing k^,, means a decrease in the 
fraction of photons last scattered from the re- ionised epochs, and so a decrease in Qp. But 
at the same time since is decreased, the electron density perturbations at last scattering 
are larger than the which partially compensates by increasing Qp. If the power 

spectrum can be approximated locally as a power law A'^{k) oc then Q^{k) oc A;^+". 

Recall that on galactic scales with k k^, n —2, while for k << k^ ~ — 1 and for 
k >> kc, n — >■ —3. So the polarisation anisotropy Q^{k) oc k'^'^"', will increse with k at 
k << kc, and decrease for >> kc, attaining a maximum at /c ~ kc- This was one of the 
reasons for normalising our estimate to k = kc- 

For small angular scales, one can also set up an approximate correspondence between 
the wavenumber k and the angular mutipole number /, using / = kR^,. Here R^, translates 
co-moving distance at the surface of last scattering, (roughly the epoch when r = r*), to 
angle on the sky, and for a flat model is given by R^, = tq — r*. In case a/ro << 1, 
this approximation will be reasonable, but for a thick last scattering surface the above 
correspondence is less accurate. In the model discussed above, where k^, = 1, a wavenumber 
k = kc = IhMpc, corresponds to / ~ 5394, and an angular scale ~ 1// ~ 0.64 arc minutes. 
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We briefly discuss the predicted small angular scale polarisation anisotropy Qp, for some 
variations on the standard CDM model, by using the same method of calculation as above. 
For example, increasing the baryon density to fib = 0.1, leads to a smaller redshift of last 
scattering with a given but also a smaller F. For = 1, one gets a slightly larger 



value yQ^{kc) ~ 2.9 x 10 ^^K, than in SCDM. Also if the primordial spectrum is tilted 
to n = 1.2, the best flt slope determined by COBE [|l^, (keeping all other parameters of 



SCDM same), then also one has a larger value \jQ'p{kc) ~ 4.0 x 10"^///^. 

Suppose we adopt a A+ CDM type model, as discussed for example in Ref. [Q, with 
VLq = 0.35, Qa = 0.65, h = 0.7, n = 1, fib = 0.04. Then we get ~ 63 for = 1. Assuming 
that is early enough that cr can be estimated as for the flat universe, we have 



a 



6000h-'Mpc frrn_^. ^^^^ 



fio^'(l + z,)i/2 Vr. 



So for the A+ CDM model, we get a ~ 687. bh'^Mpc. For this model /"H^o) ~ 1-24, 
and adopting a normalisation of the power spectrum as in [l^Jl^ we then get Ae(/Cc, z^^,) ~ 
3.56/(1 + z,), which leads to ~ 1.9 x 10-^fxK{{AT/T)Q{Te)/18fiK). Note that in 

this model the integrated Sachs- Wolfe effect, which makes a contribution to the present 
day quadrupole, will make little contribution at redshifts ~ z*. However, using Eq. (10) 
of ref. |T6[, we estimate that this will cause only a 10% reduction in the above value. For 
an open CDM model (OCDM), with A = 0, but all other parameters as for the above A+ 
CDM model, z^, is the same and Ae{kc) is of the same order, using the power spectrum as 
determined by [|l^]. However the quadrupole at last scattering is likely to be smaller; because 



the integrated Sachs Wolfe effect contributes a larger part of present day quadrupole. So 



the predicted \/Qp{kc) is likely to be smaller than the for the above models. 



B. Isocurvature type models 

Finally, consider as an alternative to the standard models, the isocurvature model re- 



cently discussed by Peebles [Q; where density perturbations are provided by CDM that 
is the remnant of a massive scalar fleld frozen from quantum fluctuations during inflation. 
The novel feature of such a picture, as pointed out in [0, is that the primeval CDM mass 
distribution is proportional to the square of a random Gaussian process; so prominent up- 
ward fluctuations are much larger (by factor F ~ 3), than for a Gaussian process with the 
same RMS. The merits of such a picture has been discussed by Peebles [|18|. We consider 
two representative models. Model 1 discussed in [|l^ adopts Qq = 0.3, A = 0.7, fib = 0.05, 



h = 0.7 and, a matter power spectrum, which, on the relevant small scales, can be approxi- 
mated as Al{k) = {k/kof-^"", with m = -1.8, ko = O.lhMpc'^. And for model 2, flo = 0.1, 
A = 0.9, fib = 0.05, h = 0.7 and, m = -1.4. 

Note that in these models, due to early structure formation, re-ionisation is expected 
to occur at large redshifts. The optical depth to electron scattering, measured from the 
present epoch could then rise to values larger than unity. However the possible ionisation 
history in these models is largely unexplored. In order to get a preliminary estimate of the 
anisotropics in polarisation that could be generated, we simply use Eq. ( |3.20|) (implicitely 
making the simplifying assumptions which went into its derivation). So the universe after 
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standard recombination is assumed to be largely neutral, and then re-ionised after an epoch 
r = r*. Again at r ~ r*, a quadrupole would arise from the free-streaming of the large 
scale entropy perturbation at recombination (the isocurvature effect cf. [jl9[). Further, in 
Eq. (p.3|) we assume to be the epoch with = 1, but take ~ 1 (to reflect the fact that 
little of the small angular scale anisotropy is due to conventional last scattering at around 
the re-combination epoch). We hope to return to a better treatment of these models in 
future work. 

For model 1 of Peebles [|1^], one then gets z^^ ~ 52, and from Eq. ( [4.5| ), a = 821h^^Mpc. 
Also for » 1, we have /(fi(2*)) ^ 1, while for = 0.3, A = 0.7, /~^(fio) ~ 1-3. 
Putting in all the numerical values, we then get 



.(AT/r)g(r,) / k 



0.1 



QUk) « 2.7F X [jj^ ) (4.6) 

Here we have scaled {AT/T)q{T(,) by a smaller value of lOfiK, since these isocurvature 
models predict a somewhat smaller quadrupole than SCDM models (cf. Fig 1. of ||18|). We 
have also incorporated a factor F, to remind ourselves that the non-gaussian statistics of the 
density field, may lead to F times larger prominent upward fluctuations. If we wish to convert 
k to I in this model, we again use = tq — t^, with tq ~ 2.17/ (fi^^iJo) ~ 862.2h~^Mpc 
(cf. Eq. (20) of ref. Q). Also r, ~ 1520.6/1"^ and so A; = IhMpc corresponds to / ~ 10, 342 



or an angular scale ~ 0.33 arc minutes. 

A similar analysis for model 2, gives z^ ~ 36, a ~ 1710h^^Mpc, and a larger polarisation 
signal, with 



.2„,^(AT/T)Q(r. 



0.3 



qm ^ 5.7F X — [jj^ j (4.7) 

Note that Peebles adopts a cosmological constant for convenience of analysis. If we were 
to consider open versions of these models, z^ and a are nearly the same (since z^, » 1, 
but the power in electron density perturbation at last scattering is much larger, because of 
a much larger f~^{Qo). The numerical value, in Eq. ( |4.6D and ( [4.7|) , at k = lhMpc~^ gets 



increased respectively to yQp{kc) ~ 4.6F x 10~'^^K (model 1) and \jQ^{kc) ~ 0.17F^K 
(model 2). The k dependence remains the same. In the open models, the density on scales 
oi k = kc are already going non-linear at and so the above numbers provide only a crude 
estimate. We see in general that these isocurvature models predict much larger polarisation 
signals compared to the SCDM type models. First, the RMS vaue is larger. Further, because 
of the non-gaussian statistics of the density field, one expects prominent upward fluctuations 
F ~ 3 times larger than the RMS value (cf. [@). 

We note in passing that the older versions of the baryonic isocurvature models, say 
with VLq ~ 0.2, VL}, = 0.05, h = 0.8, m ~ —0.5 cf. |21[], leads to even larger signals with 

~ 0.3fxK{{AT/T)Q{re)/l0f^K){k/kcY/^. However these models may already be ruled 
out by the fact that they result in spectral distortions, larger than the limit implied by the 
COBE observations [El. 
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V. CONCLUSIONS 



In this paper we have explored the possibihty of a Vishniac type contribution to the 
polarisation anisotropy at small angular scales. It is well known that non-linear effects can 
make significant contribution to temperature anisotropy on small angular scales, through 
the Vishniac effect, especially in re-ionised models. This arises due to the mode coupling 
of large angular scale, first-order velocity perturbations, with small angular scale electron 
density perturbations. We have considered here whether a similar effect contributes to the 
polarization anisotropy, by studying the coupling of large angular scale, first-order temper- 
ature anisotropy (quadrupole) with small angular scale electron density perturbations, in 
re-ionized models. 

We find that in cold dark matter models and its variants, the Vishniac type effect leads to 
a fairly small polarisation anisotropy, with \Jq^ ~ 0.02 — 0.04/xi^', on scales with k ~ IhMpc 
(or angular scales of arc minute or smaller). However in isocurvature type models the 



Vishniac type contribution can result in much larger signals. For the models of Ref. |18], the 



anisotropy on small anugular scales is non-gaussian, with prominent upward fluctuations of 
order 0.1 — 0.5 fiK, assuming F ~ 3. This reflects basically the fact that, the isocurvature 
type models have much more power on small scales and so produce much larger electron 
density fluctuations. We note in passing, that the suppression factor, due to the flnite 
thickness of the last scattering surface, on the small scale polarization anisotropy, is much 
milder than that obtained in ref. This is because, as mentioned earlier, the flrst order 
temperature quadrupole contribution, arising due to the free-streaming of the monopole at 
recombination, was not included in their analysis. 

It is clear that the polarisation signals on arc minute scales predicted by CDM type 
models will be difficult to detect, but those predicted by isocurvature type models will be 
much easier. If small scale polarisation anisotropy is eventually detected, it will open up the 
novel prospect of studying directly, both the quadrupole anisotropy and small scale electron 
denity fluctuations, at high redshifts. As pointed out in ref in a different context, one 
can then also reduce the cosmic variance of the quadrupole signiflcantly. In this paper we 
have made several approximations to analytically estimate the polarization anisotropy. We 
plan to return to a better numerical analysis in the near future. 
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